A study on the size of snow particles in powder snow 

avalanches 

Marie CLEMENT-RASTELLO* 
Cemagref, u.r. etna, 2 rue de la papeterie, BP 76, 
38402 Saint Martin d'Heres France 



Abstract 



1 Introduction 



In this work, we study the size of the par- 
ticles involved in a powder snow avalanche 
phenomenon. To determine these sizes, we 
study all the phenomena the particles have 
to face before arriving in the "body" of the 
avalanche. We study the boundary layer 
which is at the bottom of the avalanche. 
We determine, with the help of experimental 
data, the range of size of the particles that 
can be entrained by the avalanche. We then 
examine the possibility for these particles to 
reach the top of the boundary layer, and so 
to take part in the avalanche. Our final result 
is that the more frequent particles suspended 
in a powder snow avalanche have a size lower 
than 200/im. 



*This work was supported by the Pole Grenoblois 
sur les Risques Naturels 



Powder snow avalanches are very destructive 
and very ill-known phenomena. In this ar- 
ticle we try to study them with a "micro- 
scopic" scale point of view. This kind of ap- 
proach has been recently used on experimen- 
tal sites where people try to know not only 
the macroscopic but also the inner properties 
such as pressure, velocity, etc. of avalanches 
(see Dent and others (1998), Nishimura and 
others (1997), Qiu and others (1997)). The 
purpose of this study is to know a bit better 
the size of the snow particles that are encoun- 
tered in a powder snow avalanche. This has 
a real importance for someone who wants to 
study the interaction between particles and 
the turbulence of the flow. Indeed it is impor- 
tant to know the size of the particles that are 
involved in the avalanche so as to compare it 
with the characteristic scales of the turbulent 
flow, to determine the difference of velocity 
between air and particles. The study pre- 
sented here is a theoretical study in which we 
try to find out which sizes of particles are en- 
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countered more frequently in a powder snow 
avalanche. 

The way we proceed is the following: we 
study separately each of the mechanisms that 
are encountered by a snow particle before 
reaching the body of the avalanche. For each 
of them, we look at the range of sizes of the 
concerned particles. At the end we compare 
the different ranges found in each case so as 
to find the range that we must consider for 
a powder snow avalanche. The study is two- 
dimensional and deals neither with the front 
part of the avalanche nor with its tail but 
with its body. A pure powder snow avalanche 
is studied going its way on the snow cover. 
The different phenomena that are considered 
are: the size of the particles in the snow 
cover, the picking up of the particles from 
the snow cover by air friction and by colli- 
sions, and finally the proportion of each type 
of particles that manage to reach the top of 
the boundary layer, which is at the bottom 
of the avalanche, and so that take part to 
the body of the avalanche. For this study, 
we need experimental data. We chose to use 
the results obtained by Nishimura and oth- 
ers (1997) while making measurements in an 
effective avalanche. 

2 The particles of the 
snow cover 

A first range for the size of the particles is 
given by the particles which are present in 
the snow cover. This range of size is not a 
fixed parameter, it depends on many param- 



eters such as temperature, wind velocity, etc. 
during and since the snow fall. Mellor (1964) 
studied the size of the particles that can be 
found in snow layers. We use his results con- 
cerning recent snow covers. From these we 
conclude that the size of the more frequent 
particles is smaller than 1 — 2 mm. 

3 The pick up of the par- 
ticles from the ground 

3.1 Mechanisms 

A particle on the ground can be set into mo- 
tion either due to air friction or to the impact 
of another particle which is projected into the 
ground. 

Here, the aim is to obtain the size of the po- 
tentially mobilized snow grains. Let us first 
look at the collisions. What we show here is 
that the collisions between particles and the 
snow cover can increase the range of size with 
small radius but not really with large radius. 
Indeed let us first suppose the "best case" of 
energy transmission. One particle of radius 
o c arrives with a velocity u c . If it collides 
elastically with a particle of the snow cover 
of radius a/ it gives all its energy to this parti- 
cle. In this case the leaving particle will have 
a velocity 

ui = {—)%u c . (1) 

One can see that even in this ideal case to 
have a significant velocity the leaving parti- 
cle cannot have a size much larger than the 
incoming one. In reality it seems clear that 
the transmission of energy is not so perfect 



2 



and that we do not have all the energy trans- 
ferred from one particle to another. So we 
have 



Ui < 
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(2) f 



and our previous conclusion remains valid. 
On the other hand the collision phenomenon 
does not give any restriction on the lower 
limit of the sizes of the particles. So when 
we will compare the results of the picking-up 
by collisions and of the picking-up by air fric- 
tion, because of the collisions we will lose the 
information on the lower limit. For the upper 
limit the collisions will have no influence. Let 
us now look at the range we obtain from pick 
up by air friction. 
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Figure 1: Threshold friction velocity versus 
diameter of the particles 



suits they deduced the following law: 



3.2 Threshold velocity for par- 
ticle pick-up 

In order to determine the radius of the parti- 
cles which are entrained from the snow cover 
by the friction of the avalanche, we use re- 
sults that have been established for sediment 
transport by the wind and that are also used 
in snow transport by the wind. Here we use 
the results given by Iversen and White (1982) 
(see Fig. [Q. To collect their data they placed 
samples of material on a wind tunnel floor 
and determined optically the threshold fric- 
tion velocity (u*) for each of them. The parti- 
cles they used were of different kind and had 
different shapes. From their experimental re- 
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where d is 



the diameter of a particle, A is the ratio of 
the threshold friction velocity to the shallow 
water velocity of a density current and B is 
a Reynolds number. 

Iversen and White give a lower limit (B = 
0.03) but no upper limit for their law. The 
lower limit is outside the range of our study so 
does not limit us. Concerning a possible up- 
per limit we plot on the figure the curve which 
separates the region in which Iversen and 
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White made measurements and that where 
they did not make any. They say that for par- 
ticles larger than 40 \xm the above empirical 
law is generally within 5% of the experimen- 
tal values. For smaller particles it is more 
difficult to perform the experiments and so 
the precision is lower. This is of no impor- 
tance here because due to the collisions, in- 
formation on the lower limit is lost. As usu- 
ally done in snow transport studies we use 
this law here. One restriction would be that 
in their experiments Iversen and White had 
no adhesion between the particles, as we can 
have with snow. We suppose that for powder 
snow such bounds are not predominant. An- 
other restriction is that the experiments have 
been performed each time with one given size 
of particles, so that the different phenomena 
that can occur when there are both small and 
large particles are not taken into account. A 
further major restriction is that Iversen and 
White's results were established in clean air 
on a horizontal floor, whereas in our case the 
air is loaded with particles and the ground 
can be inclined. From a recent study by Bin- 
tanja (1998) it is inferred that the presence of 
particles modifies some characteristics of the 
boundary layer. However there are insuffi- 
cient results that could be taken into account 
in a modified boundary layer model. So we 
use a clear air velocity profile essentially for 
lack of better knowledge. 



3.3 Determination of the fric- 
tion velocity in a powder 
snow avalanche 



Now that a relation between the threshold 
friction velocity and the diameter of the po- 
tentially mobilized particles is available, the 
friction velocity in a powder snow avalanche 
must be determined. The avalanche stud- 
ied by Nishimura and others was a "mixed 
avalanche" : there was a dense part and above 
it a powder snow part. We suppose that 
the internal structure of the powder snow 
avalanche is not modified by the presence of 
the dense avalanche and use the results of the 
measurements for the pure powder snow. The 
velocity of the front part of the avalanche was 
from 40 m.s" 1 to 60 m.s~ l . At a given point, 
the static pressure was measured and com- 
pared to the static pressure measured some- 
where outside the influence of the flow. From 
this the velocity at the measuring point was 
measured and for the part of the avalanche 
behind the front a value of 10 m.s' 1 was 
found. The point where the measurements 
were taken was located less than seventy cen- 
timeters above the bottom of the powder 
snow part; this is in the boundary layer of the 
powder snow avalanche. Let us suppose that 
the boundary layer is of the order of one me- 
ter thick. Later on we will see that the exact 
value of the size is of no importance. Thus, 
if z < 25 cm (z < 0.255 with 8 — 1 m the 
thickness of the boundary layer), the value of 
the friction velocity can be directly estimated 
from the logarithmic law which governs the 
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velocity in this region : 

u = ^ln(-) (4) 
K z 

where K = 0.41 is the von Karman constant, 
z the roughness height, z = Co( ^ g - , with 
Co = 0.021 (see Owen (1964) and Rasmussen 
and Mikkelsen (1991)). There are, theoreti- 
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Figure 2: Determination of u* 

cally two possible values for u* (see Fig. [2]). 
In comparison with the values of u* found 
in snow transport by wind, with velocities 
comparable to those of avalanches, it appears 
that the lower value is the physically relevant 
value. We deduce from Fig. [2] that in this re- 
gion max(u*) ~ 1 m.s~ l . 
In the upper region (25 cm < z < 70 cm) the 
boundary layer profile relating velocity and 
height is unknown. If the measurement has 
been done in this region, then u(ze[0; 0.25]) < 
10 m.s -1 . So with the formula relating u 
and u* in the logarithmic region, one neces- 
sarily has u* < 1 m.s~ l . Thus here again 



max(u*) = 1 m.s^ 1 . 

Let us now suppose that 5 ^ 1 m. With the 
same reasoning as before one can show that if 
u has been measured in the upper region the 
deduced value of u* is lower than it would 
be if u had been measured in the logarithmic 
region. Because in the logarithmic region u* 
is decreasing with z (see Fig. |5J), for any 5, 
for z < 0.255 the value of max(u*) remains 
1 m.s' 1 . Thus the size of the boundary layer 
does not change anything in this estimate. 
Two things should be improved in the future : 
first a better determination of u* with more 
field data, and second to account for the pres- 
ence of the particles in the boundary layer in 
the determination of the velocity profile u(z). 

3.4 Size of the entrained parti- 
cles 

We are now able to determine the maximum 
range of sizes of the particles which can be 
picked up from the snow cover by air fric- 
tion : it is [5 * 10~ 6 m; 10~ 2 m) (see Fig. [Q. 
To determine these values of the size (d) of 
the particles, formulas given by Iversen and 
White were used (see section 1X21 . The upper 
value of d still lies in the domain of validity 
of their formulas but is outside the region in 
which they performed the experiments. Be- 
cause they did not put any upper limitation 
to their formula and because our range is not 
far from the range they studied experimen- 
tally we think their formula is applicable for 
the upper value of d. 

The upper limit of this interval is not changed 
by the collisions of particles with the snow 
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layer, but the lower limit can be changed by 
the impacts (see section 13. ip . Considering 
both air friction and collisions, the size of the 
mobilized particles is therefore d < 10" 



m. 



Sedimentation in 
boundary layer 



the 



We now study the sedimentation of the dif- 
ferent particles in the boundary layer. The 
purpose is to find out which particles "fall" 
sufficiently slowly to be present at the top of 
the boundary layer, and so to take part in the 
avalanche. 

4.1 Law relating volumetric 
concentration and height 

To find the governing law for the volumet- 
ric concentration in the boundary layer, we 
use the equations of mass conservation for the 
mean flow of air and the equation of the mean 
conservation of particles 







d(4m rj ) _ S(D, 



(5) 



(6) 



dxj dxj 

where Uj are the components of the mean ve- 
locity of air, Upj the components of the mean 
velocity of the particles, <f) the mean volumet- 
ric concentration, and D t the turbulent diffu- 
sivity that comes from the turbulent closure 

()' are turbulent 
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We suppose that the mean flow is steady and 
that the velocity and concentration profiles 
are independent of x. With the boundary 
condition that near the ground u z = (u z 
is the velocity perpendicular to the slope), 
we then have u z = in the whole bound- 
ary layer. Boundary layer theories imply that 
D t = Ku * z where a s is the Schmidt number. 
The Schmidt number is the ratio a s = — be- 
tween the diffusivity of snow particle and the 
eddy diffusivity of momentum. Some more 
advanced theories (Bintanja (1998)) take into 
account that the suspended particles modify 
the turbulence of the air, so that the diffusive 
coefficient D t is modified : 
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where R is the flux "snowdrift" Richardson 
number, p p is the density of particles, p is 
the density of air and g is the gravitational 
acceleration. The order of magnitude of the 
different quantities are the following: 
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50 m.s 



[$] w 10" 2 , [z\ w t 
and [a s ] ~ 0.5 — 1 (see section 14^21) . From this 
we can deduce that [1 + AR] ~ 1. So within 
the precision of our study, this modification 
need not to be taken into account. 
To solve the second differential equation, we 
suppose that the particle volumetric concen- 
tration and its gradient vanish infinitely far 
from the boundary layer 
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The steadiness of the flow, the independence 
of the concentration profiles of x and the fact 
that u z = imply, using equation (jSJ) 



d(-$v r cos(a) + Af 
dz 







(9) 



with a the angle of the slope and v r 
\Tt p — ~Tt\. From this and with the boui 
ary conditions (see equation E]) we deduce 
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<f) 1 is the volumetric concentration at z — z%. 
Let z\ be the roughness height, so that <f) 1 
is the volumetric concentration of the parti- 
cles near the ground, i.e., the volumetric con- 
centration of the mobilized particles. v r is a 
function of the size of the particles, so the 
variation in the boundary layer of the volu- 
metric concentration will not be the same for 
different types of particles. 

4.2 Presence of the particles 

To know which particles effectively reach the 
top of the boundary layer, we study the ratio 
between the number of particles at the top of 
the boundary layer and the number of them 
near the ground; to this end we look at the 
repartition with the radius a defined as 



R(a) = f = R 



It defines the concentration of particles by 
which the body of the avalanche is fed from 
the boundary layer. To find it the relative 



velocity v r must be determined function 
of a. We assimilate it with the sedimentation 
velocity. The sedimentation velocity is ob- 
tained from the Lagrangian equilibrium state 
where the gravity force is balanced by the 
drag force. From this equilibrium we obtain 



C d Re 2 p 
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(12) 



in which Cd is the drag coefficient and Re p = 
2ap ^ r is the particle Reynolds number. For 
spherical particles, as we suppose here the 
snow particles in the flow to be, the drag co- 
efficient can be given as follows, see Nino and 
Garcia (1994) 
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(13) 

Another parameter is the Schmidt num- 
ber (cr s ) which takes different values in 
different configurations. Householder and 
Goldschmidt (1969) found that for particles 
which are denser than the ambient fluid, the 
Schmidt number is smaller than 1. Naaim 
and Martinez (1995) performed wind tunnel 
experiments with PVC particles and found a 
Schmidt number of 0.5 to 0.6. We performed 
calculations both for a s = 1 and a s = 0.5. 
The last parameter to fix is u* . As before, 
because we look for the larger range of parti- 
cle sizes, we take u* = 1 m.s^ 1 . The results 
are reported in Fig [31 To demonstrate that 
the dependences on a s and on the thickness 
of the boundary layer are weak we plotted the 
curves for either o~ s = 1 and cr s = 0.5, and ei- 
ther S = 50 cm and S = 5 m. In all cases the 
particles which, in an important number, are 



7 




Figure 3: Repartition for a s — 1, a s — 0.5, 
5 = 50 cm and 5 = 5 m 



still present at the top of the boundary layer, 
have a size smaller than 2 * 10~ 4 m. 

5 Conclusion 

The purpose of our study was to estimate 
the size of the particles which are present in 
a powder snow avalanche. They correspond 
to the particles that reach the top of the 
boundary layer, between the ground and the 
body of the avalanche. Indeed in this way 
they can take part in the processes in the 
body of the avalanche. 

We analyzed the different phenomena which 
occur in this boundary layer and their 
influence on the size of the particles. Thus, 
we studied consecutively, the sizes of the 
particles of the snow cover, the pick-up of 
particles by the avalanche and the selection 
which is done between the particles while 
going up in the boundary layer. With this 
analyze, we were able to determine the 
maximum size of the particles which take 
part in the dynamics of a powder snow 
avalanche. These particles are those which 
have a size smaller than 200/im. This knowl- 
edge will allow to use the results, that have 
been established in diphasic flow studies, 
about the interaction between particles and 
turbulence. It will also help to improve 
numerical models which use as a parameter 
the size of the particles of snow. Finally for 
laboratory experiments, this information is 
necessary for the choice of particles. 
Another step, now, is to capture snow 
particles during the flow of a powder snow 
avalanche, and measure their size. This will 



8 



help to improve our understanding of the 
processes and our theoretical model. 
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